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Critical comments on the entropic lattice Boltzmann equation (ELBE), by Li-Shi Luo, Wei Liao, 
Xingwang Chen, Yan Peng and Wei Zhang in Ref. [3|, are based on simulations which make use of 
a model that, despite being called ELBE by the authors, is in fact fully equivalent to the standard 
lattice Bhatnagar-Gross-Krook equation. As a result, the conclusion of [l| on ELBE is circular, 
hence devoid of scientific bearing. 
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In a recent paper, authors Li- Shi Luo, Wei Liao, Xing- 
wang Chen, Yan Peng and Wei Zhang, Phys. Rev. E 83, 
056710 (2011), claim that that the entropic lattice Boltz- 
mann equation (ELBE) "does not improve the numer- 
ical stability of SRT (single-relaxation-time) or Lattice 
Bhatnagar-Gross-Krook (LBGK) model" . They go on by 
quite categorically stating that "the ELBE scheme is the 
most inferior among the LB models tested in the study, 
thus is unfit for carrying out numerical simulations in 
practice". 

In this Comment, we point out that the above state- 
ments by [l| do not bear scientific relevance. The reason 
is simple: what Luo et al implemented as ELBE, is not 
ELBE, but LBGK in disguise. 

Since the correct description of ELBE was not pre- 
sented in [l[, let us remind that in the ELBE scheme, 
populations associated with the discrete velocities Vi 
evolve according to the following kinetic equation, 

f.,ix + v„ t+l)- f,ix, t) = a/3(/^^ - /,), (1) 

In the above, f'f'^ is the local equilibrium, which mini- 
mizes the entropy function H = "^^ fi\'n-{fi/Wi), where 
weights Wi are lattice-specific. In Eq. ([1]), a is the max- 
imal over-relaxation parameter, which is operationally 
available as the positive root of the entropy condition 
(absence of subscripts denotes the full set of discrete pop- 
ulations): 



iJ(ar^ + (1 - a)/) - iJ(/) = 0. 



(2) 



This entropy estimate is key, as it assures the discrete- 
time 7?-theorem: For /3 6 [0, 1], the total entropy H{t) = 



Y,x H{f{x,t)) is not increasing, H[t -f 1) < H[t). Note 
that the validity of the if-theorem requires not just the 
equilibrium to be evaluated through the minimization of 
H but also, and most importantly, the fulfillment of the 
entropy condition ([2]). Finally, whenever the simulation is 
fully resolved (populations stay close to the local equilib- 
rium) , the maximal over-relaxation parameter a becomes 
fixed automatically to the value a = 2 [2| , and the ELBE 
dD), (12) self-consistently becomes equivalent to the LBGK 
equation and recovers Navier-Stokes equations with the 



kinematic viscosity v = c'^ 



\2p 2) 



where Cg is speed of 



sound (a 0(1) lattice-dependent constant). 

Instead of comparing with ELBE ([T]) , which, by defini- 
tion, implements a self-adjusted relaxation time through 
the entropy estimate ([2]), calculations by Luo et al make 
use of a constant relaxation time r, that is 



Ii{x 



v,,t+l)-f^{x,t) = -ifP-f,). 

T 



(3) 
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The only remaining input from ELBE in the above, is 
the local equilibrium f'['^, given by equation (14) in [l|. 
At this point, however, Eq. ([3|) is no longer ELBE, as 
it is in fact to all effects and purposes equivalent to the 
standard LBGK scheme. 

Indeed, the authors correctly admit that "one dif- 
ference between the ELBE and MRT-LBE is the 0{u^) 
terms in the odd-order equilibrium moments," and "... 
the difference in the even-order equilibrium moments ... 
is of the terms of the order O(u^). " 

This implies that the difference between the entropic 
equilibrium and its standard polynomial approximation, 
is on the order of the overall errors of the lattice Boltz- 
mann method for the low Mach number flows, and there- 
fore Eq. ^ is equivalent to the standard LBGK method. 
Since the iJ-theorem is no longer valid with (jH]), the an- 



thors could have opted to replace the equilibrium by 
a polynomial approximation, to at least second order, 
thereby completing in full the identity between the stan- 
dard LBGK and ©. 

With this assessment, and fulfilling the usual low Mach 
number restrictions, which is roughly the case for their 
simulations, the authors cannot make a "comparison of 
ELBE with LBGK", as the difference between ^ and 
LBGK is of the overall order of errors of the LB method 
for low Mach-number simulations. Therefore, results de- 
scribed as "ELBE" in [l|, should have been labeled as 
"LBGK". 

Instead, the authors continue as follows: "...based on 
our experience and understanding of the LBE, it is un- 
clear theoretically how the ELBE with a constant relax- 
ation parameter t can improve the numerical stability of 
the LBGK scheme, as it has been advocated [3,4]- " 

The reader can easily verify that neither in [3] nor [4] 
(Refs. Q and Q here), there is any claim about constant 
relaxation parameter improving stability, and whenever 
stability of ELBE was discussed in these papers, the en- 
tropy estimate ^ has always been provided (equation 
(10) in Q and equation (37) in [J]). All statements about 
ELBE in [l| are simply incorrect. 

Moreover, other independent authors, who have imple- 
mented the ELBE scheme, have made explicit use of ([2]) 
048|. It is hard to understand how/why Luo et al could 
have been misreading the papers [3|, |j| , to the point of 
attributing wrong statements to them. 

It is also informing to note that LBGK and "ELBE" 
([3]) perform so similarly that Luo ct al come to the point 
of writing that "they are so similar to each other" that 
"only the results obtained by using ELBE are shown in 
Fig. 4" (Ul: page 6). This twin-behavior alone should 
have warned the authors, that the two allegedly different 
methods were basically the same. 

On the other hand, it is clear that Luo et al. do know 
what ELBE is, as revealed at the very end of the paper. 
On page 23; the authors inform the reader that: "We did 
not test the ELBE with a variable relaxation time [52- 
54], which is supposed to guarantee numerical stability, 
because it is computationally inefficient and unphysical 
with a viscosity depending on space and time; a stable but 
inaccurate, unphysical and inefficient scheme is simply 
not a viable one. " Thus, by their own admission, what 
Luo et al have implemented is not ELBE, but just LBGK 
in disguise. 

It is not the scope of this Comment to discuss in de- 
tail the effect of the variable relaxation time in the gen- 
uine ELBE scheme, leading to the effective viscosity in 
under-resolved simulations. Still, a few clarifications are 
in order. The variable viscosity in ELBE is a built-in 
sub-grid viscosity. In the resolved direct numerical sim- 
ulation, the ELBE viscosity remains constant automati- 
cally, and corresponds to a = 2, as it has been already 
mentioned. When the grid is coarsened, local instabil- 
ities due to lack of resolution typically lead to a col- 



lapse of LBGK. It is in this situation that ELBE pro- 
ceeds with the sub-grid viscosity. To this regard, we wish 
to point out that there has never been any mystery as 
to the fact that ELBE is a natural extension of LBGK 
into sub-grid simulations; with the distinctive trait that 
the stabilization mechanism is directly informed through 
the second principle (_ff-thcorcni). This is well reflected 
by the specific way ELBE mends instabilities; most of 
the time during the simulation the relaxation parameter 
remains constant everywhere, so that indeed ELBE col- 
lapses to LBGK. It is only at the onset of a local instabil- 
ity, that ELBE deploys its built-in entropic stabilization 
capability. These stabilization events maybe rare in time 
and very localized in space (so are incipient instabilities) 
but they make the whole difference. This self-adaptive 
stabilization shows an elegance and universality (no fine 
tuning of parameters) which is simply unknown to any 
other LB method. It should be stressed that ELBE is 
not intended for fully resolved direct numerical simula- 
tion (DNS) (in this regime ELBE is identical to LBGK) 
but rather to extend LBGK to much higher Reynolds 
numbers, beyond the strict DNS regime. Although the 
test case chosen in [l| (laminar flow in a two-dimensional 
lid-driven cavity) is quite standard, on the basis of this 
test alone, it is difficult to judge any method for the pur- 
pose of DNS. Quantitative information on the adaptive 
feature of ELBE in high Reynolds number flow simula- 
tions, can be found in [5|, |9|. 

The authors of [l[ should know that most, if not all, 
under-resolved simulations, and sometimes even slightly 
over-resolved ones, imply some form of effective viscos- 
ity, which is " dependent on space and time" (the popular 
flux-corrected-transport techniques and spectral hyper- 
viscosity being two examples in point, respectively). A 
bold and categoric rejection of a "variable viscosity" as 
"unphysical", just reveals a very ad- hoc view of compu- 
tational fluid dynamics. 

As to computational inefficiency, again, there has never 
been any mystery on the fact that the implementation 
of the entropic estimate requires the solution of a single 
non-linear equation for the running parameter a, at each 
lattice site and time-step. The corresponding computa- 
tional overhead, however, is more than compensated by 
the much reduced grid demand at high Reynolds numbers 

i. 

Summarizing, ELBE with a constant relaxation time 
is not ELBE, but basically LBGK in disguise. Thus, 
what Luo et al achieve in the end, is a circular result, 
i.e. cross-compare minor LBGK variants. In the face of 
the lack of supporting evidence, in no less than 25 PRE 
pages, categoric and over-restrictive statements such as 
"the ELBE scheme is the most inferior among the LB 
models tested in the study... " appear to be totally unjus- 
tified and devoid of scientific relevance. Ref. [l| is simply 
a comparison of MRT, TRT and SRT models in a stan- 
dard flow (similar to, e. g. jlG|). 
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